Non-quantum Entanglement Resolves a Fundamental Issue in Polarization Optics 
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The issue raised in this Letter is classical in the sense of being quite ancient: which subset of 4 x 4 
real matrices should be accepted as physical Mueller matrices in polarization optics? Non-quantum 
entanglement between the polarization and spatial degrees of freedom of an electromagnetic beam 
is shown to provide the physical basis to resolve this issue in a definitive manner. 
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Entanglement is traditionally studied almost exclu- 
sively in the context of quantum systems. However, this 
notion is basically kinematic, and so is bound to present 
itself whenever and wherever the state space of interest 
is the tensor product of two (or more) vector spaces. Po- 
larization optics of paraxial electromagnetic beams hap- 
pens to have precisely this kind of a setting, and so one 
should expect entanglement to play a significant role in 
this situation. It turns out that entanglement in this 
non-quantum setup is not just a matter of academic cu- 
riosity: we shall show in this paper that consideration of 
this non-quantum entanglement resolves a fundamental 
issue in classical polarization optics. It will appear that 
this issue could not have been resolved without explicit 
consideration of entanglement. We begin by outlining 
the structure of classical polarization optics [J, 0, 0, 0] 
The Mueller-Stokes Formalism: Traditional 
Mueller-Stokes formalism applies to plane electromag- 
netic waves or, more generally, to elementary beams 
(see below). If the wave propagates along the positive 
z-axis, the components Ei, E2 of the transverse electric 
field along the x and y directions can be arranged into a 
column vector 



E \ 



El 
E2 



(1) 



called the Jones vector, analogous to the state vector of 
a qubit. [A scalar factor of the form ^^i'^^-'^t) has been 
suppressed.] While E'^ E = \Ei\'^ + \E2\'^ is {a measure 



components specifies the state of polarization. 

When E is not deterministic, the state of polarization 
is described by the coherency or polarization matrix 



(EE^) 



(EiEl) {EiE*2) 
{E2EI) (E^E^) 



(2) 



where ( ■ • • ) denotes ensemble average. The two defining 
properties of the coherency matrix are hermiticity, (f>^ — 
$, and positivity, $ > 0: every 2x2 matrix obeying 
these two conditions is a valid coherency matrix. 

It is clear that the intensity corresponds to tr<I>, and 
fully polarized (pure) states describable by Jones vectors 
E correspond to det $ = 0. Partially polarized or mixed 
states correspond to det $ > . Thus coherency matrices 
are analogous to the density operators of a qubit. 

Since is hermitian, it can be conveniently described 
as real linear combination of the four hermitian matri- 
ces To = l2x2, n = 0-3, T2 = cTi, T3 = 02 which are 
mutually orthogonal, tr TaTb ~ 25ab ■ 



1 ^ 

<i>^-^SaTa ^ Sa= tr(ra$) . 



(3) 



The coefficients Sa define the components of the Stokes 
vector S £ R^. The intensity equals So = tr $. 

While hermiticity of $ is equivalent to reality of the 
Stokes vector S G i?^, the positivity conditions tr$ > 
0, det $ > read, respectively, 6*0 > 0, Sq ~ Sf - S2 - 



of) the intensity, the ratio 7 = E1/E2 of the (complex) S§ > 0. Thus, permissible polarization states correspond 
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to the positive light cone and its interior (solid cone). 
Pure states live on the surface of this cone. 

Typical systems of interest in polarization optics are 
spatially homogeneous (in the transverse plane), in the 
sense that their action is independent of the coordinates 
(x, y). If such a system is deterministic and acts linearly 
on the field amplitude, it is described by a complex 2x2 
numerical matrix J, the Jones matrix of the system; 

J : E -> E' = JE ^ 

^ = {EE'^) ^ ^' ^ {E'E'^) ^ J^J^ . (4) 

Such Jones systems are analogous to hamiltonian evolu- 
tions of a qubit; since the intensity Sq = tr $ need not be 
preserved, J need not be unitary. It is clear that Jones 
systems map pure states (det $ = 0) into pure states. 

We can go from a pair of indices, each running over 1 
and 2, to a single index running over to 3 and vice versa. 
Thus, the elements of $ can be written as an associated 
column vector l> with $o — = *i'i2, ^2 = 'I'2i, 

and $3 = <I>22. The one-to-one relationship (3) between 
S and $ may thus be written as the vector equation 
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Thus, while $ = {EE'^) the associated column vector 
$ = {E ® E*). The 4x4 numerical matrix A exhibited 
above is essentially unitary: = \A'^ . 

Optical systems of interest can be more general than 
the ones described by Jones matrices. Such a general sys- 
tem is said to be non-deterministic, and acts directly on 
the Stokes vector rather than through the Jones vector. 
It is specified by a 4 x 4 real matrix called the Mueller 
matrix, transforming the Stokes vectors linearly: 



M : S^S'^MS. 



(6) 



Mueller matrix of a Jones system J will be called Mueller- 
Jones matrix M{J). 

Since M produces a linear transformation on S, the 
linear invertible relationship (3) or (5) between S and <& 
implies that M will induce a linear transformation 
on $, which we may write in the form : 



(M) 



(7) 



k£ 



The fact that $ ' needs to be hermitian for all hermitian 
demands that the map or super-operator If , viewed 
as a 4 X 4 matrix with ik (going over to 3) labeling the 
rows and ji labeling the rows, be hermitian. Let us define 
a new matrix J5(*^) by permuting the indices of H'^^^: 



B 



(A/) 



= H 



(M) 
ik,jl 



(8) 



i.e., _B(*^) is obtained from _H"(^) by simply interchanging 



H^i'^ with \ 

j{M) 



H, 



(M) 



03 



with H 



(M) ^(M) 
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22 



with H. 



(M) 
30 ' 



and H^^' with H^^^K In terms of B^'^''\ Eq. (7) tran- 
scribes into the vector equation $ ' = 5'*^)$, and in view 
of (5) we have the invertible relationship 



^(M) ^ A-^MA, M = 



(9) 



Thus, the correspondence between real matrices M and 
hermitian matrices H^^^\ established through (8) and 
(9), is indeed one-to-one. Elements of in terms of 

those of M can be found in Eq. (8) of Ref. [1]. 

If the system described by M is a Jones system with 
Jones matrix J, it is clear from the transformation law 
$ ^ $' = J$J^ given in (4) that B^^) = J (g, J* and, 
consequently, H'^^'^^ — Jj\ where J is the column vector 
associated with the 2x2 matrix J . Thus we arrive at 
the following result of fundamental importance |3] ■ 
Proposition 1 : A Mueller matrix M represents a Jones 
system if and only if the associated hermitian matrix 
i/*^*^-' is a one-dimensional projection. If H^'^'^ is such 
a projection Jj\ then M = M{J), J being the 2x2 
matrix associated with the column vector J. 

As a consequence, which is mathematically trivial but 
quite important for the issue on hand, we have Q 
Proposition 2 : A real matrix M can be realized as 
a positive sum (ensemble) of Mueller- Jones matrices if 
and only if the associated hermitian matrix is 
positive semidefinite. If = 'Zk ■^^''^ ^^''''^ ^ then 

M = J2k M(J('=)) where M(jW) = A{J^'''> (g> jW *)A-i 
is the Mueller- Jones matrix associated with J^''\ 

With this brief outline, we are now ready to describe 
the fundamental issue being addressed in this Letter. 
The Issue : The Mueller-Stokes formalism takes as state 
space f2(P°') the collection of all Stokes vectors: 



f^(poi) = {S eR'^ \ So>0, S^GS > } , 



(10) 



where G = diag(l, —1, —1, —1) . Thus given a 4 x 4 
real matrix M, in order that it qualifies to be a Mueller 
matrix one should demand that it maps the state space 
j^Cpoi) into itself. Let us denote by M the collection of all 
(real ) M matrices which map rj(p°i) mto itself. While the 
set of M matrices which map fi^P"'^ onto itself is clearly 
the six-parameter family SO {3, 1) U GSO{3, 1), where 
SO{3, 1) is the proper orthochronous Lorentz group, M 
is much larger; it is a sixteen-parameter family. 

Let us denote by A^^+' the collection of M matrices 
which can be realized as positive sum of Mueller- Jones 
systems M{J). It is clear that A^^"*"^ is contained in ^A. 
The structure of 7Vi(+) is fairly simple: it is clear from 
Proposition 1 fs*] that elements of Al*^"'"^ are in one-to- 
one correspondence with nonnegative 4x4 matrices [3] . 
But the structure of M is considerably more involved. 
Owing to a sequence of developments B 01 , which 
are surprisingly recent in relative terms, we now have a 
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complete characterization of ^A. The basic tool has been 
orbits of M under double-coseting by 50(3, 1) 0. 

That elements of M.'''^^ are Mueller matrices is clear, 
for they are realized as convex sums of Jones systems. 
That M matrices which fall outside M. are not Mueller 
matrices is also clear, for they fail to map the state space 
Q(poi) into itself. Thus the issue is really one about the 
grey domain Hn between', the complement of M.'^^^ in Ai: 
are these M matrices physical Mueller matrices? 

By definition, members of this domain cannot be re- 
alized as positive sums of Jones systems; but they map 
f](poi) into itself. No one has come up with a scheme 
to realize them physically. On the other hand there 
are Mueller matrices, extracted from actual experiments, 
which fall deep into this grey domain (see Ref. [1| for ex- 
amples from Ref. fl^). 

There are two difficulties in simply dismissing these 
matrices as unphysical: first, the experimenters did not 
realize them as convex sums of Jones systems, and so the 
fact that they fall outside M. cannot be enough reason 
to dismiss them; and secondly, within the Mueller-Stokes 
formalism there seems to exist no additional qualification 
we can demand of a Mueller matrix, over and above the 
requirement that it should map 17(P°1) into itself 

As a simple illustration of this grey region between 
A^(+) and Ai, let us consider Af matrices of the diagonal 
form diag(l, di, d2, ds). It is clear that M will map 
ri(P°') into f](P°'), and hence be in M, if and only if dk < 
1, /c = 1, 2, 3. In the Euclidean space spanned by 
the parameters (di, d2, d^) this corresponds to the solid 
cube with vertices at (±1, ±1, ±1). 

The hermitian matrix H^^^^ associated with 
diag(l, di, c?2, dz) is, from (10), 



1 



l + di 

l-di 

d2 - ^3 

^2 + ^3 



d2 + dz 

d2-d3 

1 - rfi 

1 + di 



(11) 



Clearly, H^'^'^ > if and only if -di - ^2 - < 1, -di + 
d2 + d^ < I, di + d2 — d^ < I, and di — d2 + d^ < 1. i.e., 
iff (di, ^2, d^) is in the solid tetrahedron with vertices at 
(1, 1, 1), (1, -1, -1), (-1, 1, -1) and (-1, -1, 1). 
Proposition 3 : For M matrices of the restricted form 
diag (1, di, c?2, ^3), Ai corresponds to the cubical region 
with vertices at (c?i, c?2, ^3) = (±1, ±1, il)j whereas 
A^(+) corresponds to the inscribed solid tetrahedron 
with vertices at (1, 1, 1), (1, —1, —1), (—1, 1, —1) and 
(-1,-1,1). 

In this Letter we present a compelling physical ground 
which judges every M matrix which is not an element 
of A4^~^'> as unphysical. And this physical ground comes 
from consideration of entanglement between the polariza- 
tion and spatial degrees of freedom of an electromagnetic 
beam. 

Non-quantum Entanglement : Let us now go be- 
yond plane waves and consider paraxial electromagnetic 



beams. The simplest beam field has, in a transverse plane 
z = constant described by coordinates {x, y) = p, the 
form E{p) — {Eix + E2y) i^{p), where Ei, E2 are com- 
plex constants, and the scalar- valued function 'ip{p) may 
be assumed to be square-integrable over the transverse 
plane: 'ip{p) £ L?{B?). x, y are unit vectors along 
the X, y axes. It is clear that the polarization part 
{Eix + E2y) and the spatial dependence or modulation 
part ?/'(p) of such a beam are well separated, allowing 
one to focus attention on one aspect at a time. When 
one is interested in only the modulation aspect, the part 
{Eix -f E2y) may be suppressed, thus leading to 'scalar 
optics'. On the other hand, if the spatial part '^(p) is 
suppressed we are led to the traditional polarization op- 
tics or Mueller-Stokes formalism for plane waves. 

Beams whose polarization and spatial modulation sep- 
arate in the above manner will be called elementary 
beams. Suppose we superpose or add two such elementary 
beam fields {ax + by) 4'{p) and {cx + dy) x{p)- The re- 
sult is not of the elementary form {ex + fy) (j){p), for any 
e, /, (j){p), unless either (a, b) is proportional to (c, d) 
so that one gets committed to a common polarization, or 
ipip) and xip) are proportional so that one gets commit- 
ted to a fixed spatial mode. Thus, the set of elementary 
fields is not closed under superposition. 

Since superposition principle is essential for optics, we 
are led to consider beam fields of the more general form 
E{p) = Ei{p)x+E2{p)y, and consequently to pay atten- 
tion to the implications of inseparability or entanglement 
of polarization and spatial variation. This more general 
form is obviously closed under superposition. We may 
write E(p) as a (generalised) Jones vector 



El 



Eiip) 
E2{p) 



E,{p), E2{p)€L\R^). (12) 



The intensity at location p corresponds to |i?i(p)p -I- 
\E2{p)\'^ . This field is of the elementary or separable 
form iff Ei{p) and i?2(p) arc linearly dependent (pro- 
portional to one another). Otherwise, polarization and 
spatial modulation are inseparably entangled. 

The point is that the set of possible beam fields in 
a transverse plane constitutes the tensor product space 

® L^(i?^). But the set of all elementary fields consti- 
tutes just the set product x L'^{B?) of and L'^{B?), 
and hence forms a measure zero subset of the tensor prod- 
uct (g) L^{R'^). In other words, in beam fields repre- 
sented by generic elements of ® L^{R?) polarization 
and spatial modulation should be expected to be entan- 
gled : Entanglement is not an exception; it is the rule 
in ® L'^{R?), the space of pure states appropriate for 
electromagnetic beams. 

To handle fluctuating beams, we need the beam- 
coherence-polarization (BCP) matrix $(/>; p') = 
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{E{p)E{py)^: 



{E,ip)E,{p'r) {E,{p)E2{p') 
{E2{p)E^{p'T) {E2{p)E2{p') 



(13) 



As the name suggests, the BCP matrix describes both 
the coherence and polarization properties. It is a gen- 
erahzation of the numerical coherency matrix of Eq. (2), 
now to the case of beam fields. 

It is clear from the very definition (13) of BCP matrix 
that this matrix kernel, viewed as an operator from ® 
L'^{R^) ^ (g) L^(i2^), is hermitian nonnegative: 

$,fe(p;p') = 'i>fc,(p';p)*, J, fc = 1, 2; 
j d'pd'p'E{py^p; p')E{p') > 0, (14) 

VE{p) e L'^{R^). These are the defining proper- 
ties of the BCP matrix: every 2x2 matrix of two-point 
functions ^jk{p'-, p') meeting just these two conditions is 
a valid BCP matrix of some beam of light. 
Resolution of the Issue: In the BCP matrix (16), 
each of the four blocks <&jfc(p; p') = {Ej{p)Ek{p')*) is an 
(infinite-dimensional) operator L^(R^) L^(i2^). For 
the issue on hand, however, it turns out to be sufficient to 
restrict the spatial dependence to just two orthonormal 
spatial modes ?Ai(p), ip2ip)- This amounts to considering 
in place of L^{R'^) the two-dimensional space C^, the 
linear span of ipi{p) and ip2{p), so that the BCP matrix 
^(p; p') is a positive operator mapping C^(8)C^ C'^®C'^, 
the first being for the polarization degree of freedom 
and the second for the spatial degree of freedom. 

Choice of a product basis in ® transcribes the 
BCP matrix into a hermitian nonnegative numerical 4x4 
matrix of four 2x2 blocks, with the polarization (Roman) 
indices labeling the blocks and the spatial (Greek) indices 
labeling entries within each block 1^. The following or- 
thonormal product basis suggests itself naturally: 



Xia(p) 



V'a(p), X2a(p) 



Mp), (15) 



for a = 1, 2. Thus a Jones vector E{p) in (8) nec- 
essarily has the form J2j=o X]a=o ^jaXjaip)- It is of the 
separable or elementary form if and only if Cn/C 12 — 
C21/C22. It can be identified with a four-dimensional 
numerical column vector C whose entries are the ex- 
pansion coefficients Cja- C ~ [Cn, C12, C21, C22]"^. It 
follows that the corresponding (pure state) BCP matrix 
^{p,p') = E{p)E{p') can be identified with the (numer- 
ical) projection matrix $ = CC^. The M matrix acts 
on $ through the associated hermitian matrix : 



H 



(M) 



(M) 



(16) 



Note that the Greek indices are left unaffected, consistent 
with the (transverse) spatial homogeneity of M . 



Now consider the special Jones vector E{p) = Xii(p) + 
X22(p) corresponding to C = [1, 0, 0, 1]-^ and $ = 
CC^. Since Cja = Sja-, we have ^ja,kp = ^ja^kp- That 
is, the only nonzero elements of the 4x4 matrix $ are 
$00 = $03 = $30 = $33 = 1- 

Given an M matrix, let be the result of the action 
of M on $. We have from (16) 



= H. 



[M) 
ia.JP 



(17) 
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This is our final result. Suppose is not positive then 

the result of M acting on the entangled Jones vector 
E{p) = Xii(p) +X22(p), fails to be positive and hence is 
unphysical, showing in turn that M could not have been 
physical. Thus H'-^^^ > is a necessary condition for M 
to be a physical Mueller matrix. On the other hand, we 
have seen that if > then M can be realized as a 

convex sum of Jones systems, showing that _ff > is 
a sufficient condition for M to be a Mueller matrix. We 
thus have 

Theorem : The necessary and sufficient condition for AI 
to be a physical Mueller matrix is that the associated her- 
mitian matrix > 0. Every physical Mueller matrix 
is a convex sum of Mueller- Jones matrices. 

Thus M matrices which map r2(P°i) into itself should 
be called pre-Mueller matrices rather than Mueller matri- 
ces. For, to be promoted to the status of Mueller matrices 
they needs to meet the stronger condition > aris- 

ing from consideration of entanglement. In view of the 
rapidly growing current interest in an unified approach to 
coherence and polarization in optics, it is hoped that our 
result will stimulate further research into the kinematic 
role of entanglement in classical optics. 
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